Influence of O2 and N2 on the conductivity of carbon nanotube networks 
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We have performed experiments on single-wall carbon nanotube (SWNT) networks and compared with 
density-functional theory (DFT) calculations to identify the microscopic origin of the observed sensitivity of 
the network conductivity to physisorbed O2 and N2. Previous DFT calculations of the transmission function for 
isolated pristine SWNTs have found physisorbed molecules have little influence on their conductivity However, 
by calculating the four-terminal transmission function of crossed SWNT junctions, we show that physisorbed 
O2 and N2 do affect the junction's conductance. This may be understood as an increase in tunneling probability 
due to hopping via molecular orbitals. We find the effect is substantially larger for O2 than for N2, and for 
semiconducting rather than metallic SWNTs junctions, in agreement with experiment. 

PACS numbers: 73.63.-b, 68.43.-h, 73.50.Lw 



I. INTRODUCTION 

Using single-wall carbon nanotubes (SWNTs) as nanosen- 
sors, both individually and in SWNT networks, has been 
one of the most promising potential applications of SWNTs 
since their discovery. ^'^ Several experimental studies have 
demonstrated that the conductance of SWNT systems is rather 
sensitive to the presence of even single-molecule concentra- 
tions of physisorbed gas molecules such as O2 and N2.^~^ 
Further, by measuring conductivity of individually charac- 
terized SWNTs, as well as thick (metal-like) and thin 
(semiconductor- like) SWNT networks, ^'^'^^ the response of 
SWNTs to contaminants has been shown to correlate with the 
intrinsic electronic properties of the material. For example, 
it has been found that the presence of I0W-O2 concentrations, 
independent of temperature, introduces an increase in conduc- 
tance of approximately 20% on thin SWNT networks, while 
an increase in conductance of only about 1% is found for thick 
SWNT networks.^ 

On the other hand, previous theoretical studies have found 
that SWNTs are rather inert, so that gases tend only to ph- 
ysisorb to the SWNT surface. For this reason, it was 
suggested that O2 should not effect conductance through 
SWNTs, but only influence conductance at either SWNT- 
SWNT junctions, at the SWNT-metal contacts, or at SWNT 
defect sites. ^^'^^ Although the conductivity of SWNTs with 
molecules physisorbed at defect sites has been extensively 
studied,^^-^^ the conductivity of four-terminal SWNT-SWNT 
junctions has been previously studied only for small pristine 
metallic SWNTs. ^^'^^ The possible influence of physisorbed 
molecules on SWNT-SWNT junctions has not been investi- 
gated. 

In this paper we address the microscopic origin of the in- 
crease in conductance of SWNT networks when exposed to 
O2 or N2 gas. To this end, we have performed density- 
functional theory (DFT) calculations of the intratube trans- 
mission within a SWNT and the intertube transmission be- 
tween two SWNTs in the nonequilibrium Green's function 
(NEGF) formalism for O2 and N2 molecules physisorbed in 



(7,7) metallic armchair, (12,0) semimetallic zigzag, and (13,0) 
semiconducting zigzag SWNT junctions, shown schemati- 
cally in Fig. 1. Comparing our theoretical results for SWNT 
junctions with experimental measurements for SWNT net- 
works suggests that the surprising sensitivity to O2 and N2 
may be partially due to an increased tunneling probability 
through O2 and N2 physisorbed at SWNT junctions. 

In Sec. II we describe experimental measurements of the 
influence of both O2 and N2 on the conductivity of SWNT 
networks and the characterization of these networks using Ra- 
man spectroscopy. A description of the DFT and NEGF model 
used to describe the microscopic origin of this effect is then 
provided in Sec. III. In Sec. IV we compare our theoretical 
results for the SWNT junction transmission with the SWNT 
network experiments, followed by a concluding section. 



II. EXPERIMENTAL RESULTS 

Below we give a brief discussion of our experiments on 
SWNT network conductivity. A more detailed description 
may be found in Refs. 8 and 9. Figure 2 shows experimen- 
tal measurements of the conductance sensitivity to O2 and N2 
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FIG. 1: (Color online) Schematics of a (13,0) SWNT junction with 
(a) physisorbed O2 and (b) physisorbed N2. 
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FIG. 2: (Color online) Fractional change in conductance AG = 
G/G{0) — 1 versus time t in seconds and hours (inset) following 
exposure to O2 and N2 for thin (□, ) and thick (O, ) SWNT net- 
works respectively, on log-log and linear (inset) scales (Refs. 8 and 
9). 



exposure for thick and thin SWNT networks. These samples 
were initially placed under vacuum (~lxlO~^ mbar) and ir- 
radiated by a UV light-emitting diode (LED) (A ^ 400 nm) at 
low intensity (~0.03 mW/cm^) for approximately 12 h to des- 
orb surface and interbundle adsorbates (surface dopants) from 
the SWNTs. Once the SWNT network's conductance stabi- 
lized, the samples were exposed to either O2 (99.5% pure) or 
N2 (99.998% pure) at 1 atm. The conductance of the sam- 
ples was then monitored by periodically sampling (At ^ I s) 
the current while applying a fixed bias of 1 mV to the thick 
(metal-like) SWNT network (R ^ 1 kQ) and 10 mV to the 
thin (semiconductor-like) SWNT network (R ^ 1000 kQ), as 
shown in Fig. 2. 

After 5 min of exposure to O2, the thin network shows an 
increase in conductance of about 13% while the thick net- 
work's conductance changes by about 7%. For the same ex- 
posure to N2, both networks show substantially smaller con- 
ductance changes of 2%-3%. However, at exposure times of 
more than 2 h, the thin SWNT network response to N2 is sim- 
ilar to that of the thick SWNT network to O2. This might be 
caused by a weaker physisorption of N2 to the SWNT net- 
works than O2. The inset of Fig. 2 also shows that at very 
long exposure times the fractional change in conductance, 
AG = G/G(0) - 1, becomes saturated after 24 h. Further, the 
response to O2 depicted in Fig. 2 shows that the conductance 
change for a thin SWNT network is about two to three times 
that of the thick SWNT network at all times. This suggests 
that the conductance change under O2 exposure is an intrinsic 
property of the SWNT networks, present even at very low O2 
concentrations. Herein we shall focus on the microscopic ori- 
gin of the network sensitivity to O2 and N2, with the temporal 
behavior of the networks discussed elsewhere. ^'^ 

We have performed the Raman spectroscopy to character- 
ize our SWNT network samples, which were produced via 
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FIG. 3: (Color online) Raman spectra and approximate diameter dis- 
tribution of HiPco SWNT sample for an excitation wavelength (top) 
Aexc ~ 532 nm (lower black curve), Aexc ~ 785 nm (upper red curve), 
and (bottom) Aexc ~ 632.5 nm ( ). The DFT calculated diame- 
ters of d ^ 9.76 A, 9.79 A, and 10.66 A for (7,7), (12,0), and (13,0) 
SWNTs, respectively, are provided for comparison (dashed lines). 



the high-pressure carbon monoxide (HiPco) method. Fig- 
ure 3 shows the radial breating mode (RBM) Raman sig- 
nals of HiPco samples at excitation wavelengths Aexc ~ 532 
nm, Aexc ^ 632.5 nm and Aexc ^ 785 nm. The van Hove 
singularity energy separation was calculated using the tight- 
binding approximation with the carbon-carbon interaction en- 
ergy 7o 2.9 eV and carbon-carbon bond length ac-c ~ 1 .44 
A. The SWNT diameter d dependence of the RBM frequency 
^RBM for isolated SWNTs on Si02 has been shown^^ to be- 
have as i^RBM ^ 248/(it. The DFT calculated diameters for 
(7,7), (12,0), and (13,0) SWNTs ofd^ 9.76, 9.79, and 10.66 
A, respectively, are found to correlate well with the HiPco Ra- 
man shift, as shown in Fig. 3. This should ensure a good de- 
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FIG. 4: (Color online) Schematic of (13,0) SWNT junction consist- 
ing of four SWNT leads (light gray) coupled to the central scattering 
region via four SWNT primitive unit cells or layers (gray) "frozen" 
at their relative positions in the isolated SWNT. The atomic positions 
of the remaining eight SWNT layers (dark gray) and the physisorbed 
O2 molecule have been relaxed. A SWNT separation consistent with 
experiments of approximately 3.4 A has been used. The intratube 
transport through the SWNT and the intertube transport between the 
SWNTs has been shown schematically. 



scription of the SWNT network's work function, which may 
be significantly different for smaller tubes. 



III. BASIC THEORY 

Our DFT calculations have been performed with the SIESTA 
DFT code^^'^^ using a double-zeta polarized (DZP) basis set 
for the physisorbed molecules (O and N), and a single-zeta 
polarized (SZP) basis set for the SWNTs (C), and the Perdew- 
Burke-Ernzerhof exchange-correlation functional. We note 
here that the DZP and SZP SIESTA basis sets have recently 
been shown to yield transmission functions in quantitative 
agreement with plane-wave codes and maximally localized 
Wannier f unctions. When modeling O2 we have performed 
spin-polarized calculations^^ but have performed spin- unpo- 
larized calculations otherwise. 

We have modeled the junction system using 6(11) primi- 
tive unit cells or layers for each zigzag(armchair) SWNT per 
supercell, with a separation of approximately 3.4 A, as de- 
picted for a (13,0) SWNT junction in Figs. 1 and The four 
SWNT layers at the boundaries of the central region, shown 
in gray in Fig. 4, were kept fixed at their relaxed positions in 
the isolated SWNT. At the same time the central 4(9) primi- 
tive unit cells from each tube, shown in dark gray in Fig. 4, 
and the physisorbed molecules were relaxed until a maximum 
force of less than 0.1 eV/A was obtained. Since the super- 
cell has dimensions of > 25 A for each SWNT junction, a F 
point calculation was sufficient to describe the periodicity of 
the structure. 

Such a large supercell was necessary for the Hamiltonian 



TABLE I: Change in intertube conductance ACinter at the valence- 
band maximum £vb relative to the pristine junction for O2 and 
N2 physisorbed in SWNT junctions of (7,7), (12,0), and (13,0) 
SWNTs, with the respective SWNT-O2 and SWNT-N2 separations 
d in Angstrom. 
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of each of the four SWNT layers adjacent to the boundaries 
iJ^^, to be within 0. 1 eV of the Hamiltonian for the respective 
leads i^Q,, so that max \ — Ro\ < 0.1 eV. In this way 
the electronic structure at the edges of the central region was 
ensured to be converged to that in the leads. 

The Landauer-Biitticker conductance for a multi- terminal 
system can be calculated from the Green's function of the cen- 
tral region, Gc, according to the formula^^"^"^ 



out I 



(1) 



where the trace runs over all localized basis functions in the 
central region. To describe the conductance at small bias 
for semiconducting systems, the Fermi energy ep should be 
taken as the energy of the valence-band maximum syb or 
conduction-band minimum Ecb for p-type and n-type semi- 
conductors, respectively. The central region Green's function 
is calculated from 



Gc{e) 



zSc - He -Y.^a{e) 



(2) 



where z = e + i0+, Sc and He are the overlap matrix and 
Kohn-Sham Hamiltonian matrix of the central region in the 
localized basis, Sq, is the self-energy of lead a. 



Hhj. (3) 



and the coupling elements between the central region and lead 
a for the overlap and Kohn-Sham Hamiltonian are Sca and 
Hca respectively. 

The coupling strengths of the input and output leads are 

then given by Tin/out = ^(^in/out - ^In/out)- ^ f^^^- 
terminal SWNT junction, the intratube and intertube trans- 
mission functions are calculated by choosing the appropriate 
output lead, as depicted in Fig. 4. 



IV. THEORETICAL RESULTS & DISCUSSION 

For each of the three types of SWNT junctions considered, 
we find that both O2 and N2 are physisorbed with binding en- 
ergies of ~ 0.2 eV, as depicted in Fig. 1. Further, the SWNT- 
02 and SWNT-N2 equilibrium separation distance d is in the 



4 



N2HOMO 02H0M0s(T) 02LUM0sa) 




Energy (eV) 

FIG. 5: (Color online) (a)-(c) Intratube transmission and (d)-(f) intertube transmission vs energy in eV relative to the Fermi energy for a 
pristine (black solid curves) SWNT junction, and with N2 (blue dashed curves), and O2 (red dash-dotted curves) physisorbed consisting of [(a) 
and (d)] metallic (7,7) SWNTs, [(b) and (e)] semimetaUic (12,0) SWNTs, and [(c) and (f)] semiconducting (13,0) SWNTs. Schematics and 
isosurfaces of ±0.02e/A^ for molecular orbitals on the physisorbed molecules are shown above the corresponding eigenenergies (dotted lines). 



range 2.3 — 2.8 A, as given in Table I. These results agree qual- 
itatively with previous theoretical studies for O2 binding dis- 
tances and energies on isolated SWNTs. 

Figures 5(a)-5(c) show the intratube transmission for three 
prototypical SWNTs commonly found in experimental HiPco 
samples, as shown in Fig. 3. In Fig. 5(a) we see that for a 
metallic armchair (7,7) SWNT, transmission occurs through 
two channels at the Fermi level. We see in Fig. 5(b) that 
the conductance for the semimetallic zigzag (12,0) SWNT 
resembles that found in Fig. 5(a) for a metallic SWNT, ex- 
cept for a tiny band gap of < 0.05 eV at the Fermi level. In 
Fig. 5(c) we find for a semiconducting zigzag (13,0) SWNT 
a band gap of approximately 0.6 eV between the valence 
and conduction bands, through which no transmission oc- 



curs. This is only slightly smaller than the expected band 
gap of ~ 0.7 eV, based on a d"^ fit to experimental data.^^ 
These results for the intratube transmission of pristine SWNTs 
also agree qualitatively with previous DFT studies of isolated 
(5,5), (10,10), (10,5), (11,0), and (12,0) SWNTs.^^'^^'^^ We 
also find in Figs. 5(a)-5(c) that neither O2 nor N2 physisorbed 
at a SWNT junction noticeably influence the intratube trans- 
mission. 

In Fig. 5 we also show isosurfaces and eigenenergies for 
the highest occupied and lowest occupied molecular orbitals 
(HOMO and LUMO) on physisorbed O2 and N2. For these 
weakly coupled molecules, the renormalized molecular levels 
may easily be identified with the molecular orbitals of the free 
O2 and N2 molecules. Since the position of the molecular 
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levels is rather insensitive to the type of junction, it should 
also be insensitive to the exact binding geometry. This sug- 
gests that additional physisorbed molecules will influence the 
intertube transmission similarly. 

We find the intertube transmission is proportional to the 
density of states (DOS) for the system with peaks in the 
transmission at the van Hove singularities. This is consistent 
with transport between the SWNTs occurring in the tunneling 
regime, as expected for a SWNT separation of approximately 
3.4 A. 

The presence of physisorbed molecules in the SWNT- 
SWNT gap should then increase the tunneling probability at 
energies near the eigenenergies of the molecular orbitals. This 
is evidenced by the distinct peaks in the intertube transmis- 
sion for each SWNT junction at energies corresponding to the 
HOMO of N2 and the spin polarized HOMOs and LUMOs 
of O2, as seen in Figs. 5(d)-5(f). Under such conditions, a 
SWNT junction behaves as a simple tunneling electron mi- 
croscopy (TEM) tip. By applying appropriate bias voltages, 
one may potentially probe the molecular orbitals of a ph- 
ysisorbed molecule to determine its chemical composition. 

For this reason, the sensitivity of SWNT network conduc- 
tivity to O2 is at least partly due to the close proximity of the 
O2 HOMOs to the Fermi energies of typical SWNTs 0.6 
eV), as shown in Fig. 5. Further, it has been shown experi- 
mentally that defects inherent in physically realizable SWNTs 
yield p-type semiconductors. ^'"^^'"^^ The conductivity measured 
experimentally at small bias is thus at the energy of the va- 
lence band ^vb • 

As seen in Fig. 5(f), the O2 HOMO eigenenergy is only 
about 0.3 eV below ^vb for a semiconducting (13,0) SWNT 
junction. As shown in Table I, this yields a substantial in- 
crease in the intertube conductance at zero bias for semicon- 
ducting junctions in the presence of O2, while much smaller 
increases are found for the metallic and semimetallic junc- 
tions, in agreement with experiment.^ On the other hand, we 
also find physisorbed N2 increases the intertube conductance 



only slightly, also in qualitative agreement with experiment,^ 
as shown in Fig. 2. 

Although it is well-known DFT calculations underestimate 
band gaps,^^'^^'^^ since we are primarily interested in how the 
presence of O2 or N2 qualitatively changes the DOS and con- 
ductance, such calculations are still useful. 



V. CONCLUSIONS 

In conclusion, we have proposed a possible microscopic ex- 
planation for the experimentally observed sensitivity of the 
electrical conductance of carbon nanotube networks to oxy- 
gen and nitrogen gases. Our DFT calculations suggests that 
O2 and N2 physisorbed in crossed SWNT junctions can have 
a large influence on the intertube conductance. In particular, 
for O2 the close proximity of the highest occupied molecular 
orbitals with the Fermi level of the SWNT significantly in- 
creases electron tunneling across the gap. The effect is found 
to be larger for O2 than for N2 and for semiconducting rather 
than metallic SWNTs, in agreement with the experimental ob- 
servations. Our results suggests that the electrical properties 
of SWNT networks are to a large extent determined by crossed 
SWNT junctions. 
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